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ABSTRACT 

We review recent progress in the description and understanding of disoriented 
chiral condensates. Certain important unsolved issues are underlined, and the pre- 
liminary results of our program of investigation of these issues in the framework 
of the classical linear sigma model are reported. We also briefly review a formal- 
ism which could be useful at the full non-equilibrium quantum field theory level of 
analysis. 
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1 Introduction 

Recently, in order to explain rare events with a deficit or excess of neutral pions 
observed in cosmic ray experiments, there has been increased interest in the con- 
jecture that it might be possible to produce disoriented chiral condensates (DCCs), 
i.e. correlated regions wherein the quark condensate, (0\qLQR\0), is chirally rotated 
from its usual orientation in isospin space. 

On the theoretical side there has been great interest (see, e.g., Refs. [1-14]) 
both in the development of technical tools suitable for the description of this new 
phenomenon, and in the exploration of the possibilities opened by DCCs as probes 
of the structure of chromodynamics, most notably in relation to the chiral phase 
transition. 

The idea that such DCCs might be produced in high energy collisions at exist- 
ing or planned hadron or heavy-ion accelerators has originated several experimen- 
tal proposals. In particular, one of us is co-spokesman for a Fermilab experiment 
looking for DCCs in hadron-hadron collisions. In high energy p-p collisions 
which lead to a sizeable multiplicity of produced particles, but not necessarily with 
high-pr jets in the final state, the time evolution is quasi-macroscopic, because the 
hadronization time can be rather large, 3-5/. At times t before hadronization, the 
initial state partons, produced in a volume much smaller than a cubic fermi, stream 
outward at essentially the speed of light in all directions, occupying the surface of 
a sphere of radius t (in units such that the speed of light is 1). Most of the out- 
going energy /momentum is expected to be concentrated near the light cone, i.e. 
on the fireball surface. However, the interior of the fireball is also an interesting 
place. If its energy density is low enough, the interior should look very similar to 
the vacuum, with an associated non-vanishing quark condensate. Since the energy 
density from the intrinsic chiral symmetry breaking is small || [17) , and the fireball 
surface isolates the interior from the exterior of the light cone, it is reasonable to 
consider the possibility that well inside the light cone the quark condensate might 
be chirally rotated from its usual orientation. At late times this disoriented vac- 
uum would relax back to ordinary vacuum, radiating its collective modes, pions. 
The properties of the radiated pions would be strongly affected by the semiclassical, 
coherent nature of the process; in particular, one could expect anomalously large 
event-by-event fluctuations in the ratio of the number of charged pions to neutral 
pions produced. Importantly, assuming that the event-by-event deviation of the 
quark condensate from its usual orientation be random, one finds [1-4,16,17] that 
the distribution P(f) of the neutral fraction 

f = N *° = N *° 

is given by 

dP _ 1 
~¥~ V7 ' 

at large N tot . Most notably, this implies that for "DCC pions" the probability of 
finding extreme values of / is very different from ordinary pion production (which 
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is given by a binomial distribution), in which the fluctuations are expected to be 
peaked at / = 1 /3 and fall exponentially away from the peak. Recent experimental 
DCC searches are largely based on the structure of Eq. (0). 



2 Description at the Classical Level 

Since the DCC is essentially a classical pion field, the description of its space- 
time evolution is most naturally described using classical or semiclassical techniques. 
The 0(4) a- model is typically used as a model of chromodynamics in DCC studies. 
The <7-model is simple enough to be treatable, has the correct chiral symmetry 
properties, and describes the low energy phenomenology of pions. The Lagrangian 
of the (linear) cx-model is (in the chiral limit m^ — ^S) 

£ = \{d»a? + \(d^f - \{a 2 + n 2 - Jlf . (3) 

A meaningful cx-model description can start at some small proper time, of or- 
der 0.2-0.3/, near the light cone, when the collective coordinates a and n become 
relevant JCT)), [TJJ]. At this early proper time the distribution of the chiral field 

$ = CX + 27TT , (4) 



can be expected to be noisy, but with ($) =0. 

As proper time increases the field $ rolls into a minimum with $ + $ = f 2 , and 
during this "rolling phase" the pion mass is imaginary, leading to unstable growth 
of the Goldstone modes [§, 10| . Since, as mentioned in the Introduction, the energy 



density from the intrinsic chiral symmetry breaking is small || [17| , and the fireball 
surface isolates the interior from the exterior of the light cone, it is reasonable to 
expect that the interior ends up in a disoriented vacuum. 

At late times one such region of disoriented vacuum with a given isospin orien- 
tation would relax back to ordinary vacuum ($) = (ex) =f 7T , radiating pions with the 
same isospin orientation. 

In modelling these stages of evolution, the chiral limit can be safely taken as 
long as the proper time is small compared to m." 1 , while at times of order l-2m~ 1 
the pion mass can no longer be neglected and one should JK| decompose the DCC 
field into physical-pion normal modes and let them propagate out to infinity as free 
states. 

It is also reasonable to expect that approximations based on the replacement of 
the full linear cx-model @ by the simpler non-linear cx-model 

C = \{d,af+ l -{d^)\ with a 2 + vf 2 = f I , (5) 

could be reliably used at times late enough for the chiral field to have already rolled 
into a minimum with $ + $ = o" 2 + n 2 = f^. 



Actually the entire picture of DCC evolution given above can be implemented 
(although, at least at early times, the quantitative aspects are not accurately evalu- 
ated) within the non-linear cr-model in the framework of the set of classical solutions 
identified by Anselm |l] and others 0, |17| , which have the form 



$ = UV+e w ^V R , (6) 

where Vl and Vr are constant but otherwise arbitrary matrices U, and 9 is such 
that 

□0 = . (7) 

In particular, by taking Vl = Vr and using the following solution of (^) 

tt 9T 4- r — t 
6{r) = — Q(2T + r-t)e(t-r) 

jv 2T — r — t 

-Q(2T-r-t)Q(t + r) , (8) 

where B is the step function, one can describe several stages of the evolution of an 
ideal spherically symmetric DCC formed at the collision centre. For t < 2T the 
DCC field described by (^) expands at the speed of light (mimicking the "rolling 
phase" of the DCC evolution). For 2T <t< 4T the DCC field keeps expanding but 
the true vacuum starts breaking into the interior of the light cone. For t > 4T the 
true vacuum is everywhere apart from a bump of DCC field propagating outward. 
The final stage of evolution, the one in which the classical DCC field radiates pions, 
cannot be reproduced within (H); it requires, as stated above, taking into account 
the non- vanishing pion mass. 

There are two more aspects of the simple solution (§) which are worth empha- 
sizing. First we notice that (§) is a solution of (^) in the presence of a DCC source 
term on the light cone at small times; such a source term is a straightforward way 
to mimic the unstable growth of the Goldstone modes associated with the "rolling 
phase" , and could also be useful in more detailed analyses. [Note that taking Vl = Vr 
is very important, since otherwise the source term on the light cone will not vanish 
at large times.] A second aspect of (S) which could be instructive for other de- 
velopments is the dominant role played by the Minkowski geometry. This remark 
also gives us a chance to state our belief that geometry should play a central role 
in DCC physics, and therefore the investigation of the possibility of DCC forma- 
tion in hadron-hadron collisions should proceed differently from the corresponding 
heavy-ion case. 

Going back to the study of the space-time evolution of the DCC, let us now 
consider the full linear cr-model, which is necessary for a description starting at 
early proper times. A good starting point is the generalization to the linear cr-model 
of the solutions ([]) and ([F|) of the non-linear cr-model. One finds solutions of the 
form Jl3 



$ = pV + e t6TS V , 



(9) 



with 6 and p such that 



ne = o, □p- p (^) 2 + Ap(p 2 -^) = 0. (10) 

The equation for p is not easy to handle analytically, but progress can be achieved 
by numerical techniques. 

We are at present using numerical and analytic techniques in a study of the 
space-time evolution of DCCs within the linear a- model ([$]). The simulations start 
from a given configuration for the pion/sigma fields at an initial proper time, and 
we observe the dynamics of the transition to the ordered broken symmetry phase as 
the system expands. We plan to study different types of initial conditions, with the 
chiral field starting either from the top or from the brim of the "Mexican hat" and 
with various assignments for the initial velocities in the 7? and a directions. One of 
the issues that we are interested in investigating is the possibility of correlations in 
the evolution of two (or more) pieces of DCC in different regions of the lego plot. It 
would be interesting to check whether these pieces of DCC "attract" or "repel" , and 
whether the presence of DCC in one region of the lego plot stimulates production in 
neighbouring regions. Ultimately, we hope to create simulations realistic enough to 
be used in the interpretation of hadron-hadron experimental results relevant for the 
DCC, but at present our research program is still in its preliminary stages. We have 
started by investigating the linear cr-model evolution of a spherically symmetric 
DCC formed at the collision centre, and observed that, as expected W[ M, the 



evolution proceeds through the same stages of the simple non-linear a-model solution 
associated with ([8]); however, also depending on the initial conditions, we find some 
evidence of a possible secondary wave front resulting from implosion of the true 
vacuum as it occupies the interior region of the fireball. Such a secondary wave 
could in principle be very important for experimental purposes, but much more 
numerical work is needed to check whether in fact it can be expected to occur. Other 
preliminary results concern simulations of scenarios with two pieces of DCC and, 
for certain initial conditions, they indicate that the evolution of one piece of DCC 
is substantially affected by the presence of the other piece. Also for this effect much 
more numerical work is needed in order to establish its significance for experimental 
purposes. Our analytic and numerical work is also aimed at describing, within the 
proper framework of Minkowskian geometry, the way a piece of DCC affects the 
evolution of another piece of DCC by acting as a source. 



3 Quantum Field Theory 

The expression in quantum terms of the ideas discussed classically in the previous 
section is most natural in the framework of coherent states or squeezed states. We 
refer the reader to Refs. || || [l6|, 17 1 for analyses within these formalisms. Here we 
will briefly review a formalism that could be useful in investigating whether such 
quantum states are actually created as a consequence of the relevant non-equilibrium 
quantum dynamics. We are not aware of much progress in the direction of exploiting 
this formalism in the investigation of DCCs, at least in association with hadron- 



hadron collisions^, but a few conceptual points can be made even before actual 
calculations. 

One can describe a physical system via its density matrix p, p = e~^ H j (tr e~^ H ), 
and average values of observables O are determined by the density matrix: (O) = 
trpO. In equilibrium physics the density matrix is time-independent; in general, 
however, the density matrix is time-dependent, and the task of non-equilibrium 
quantum field theory is to study the time evolution of p. Non-equilibrium quantum 
physics is a vast subject, and there is no canonical approach to its investigation. 
Usually the approach is suggested by the specifics of the physical system that one 
wants to describe. We review an approach |2l], [22|] that has proved useful in early 
Universe cosmology fl23| , and might also be useful in the investigation of DCCs. This 
approach is set up in the framework of the field theoretic Schrodinger picture, which 
is particularly suitable to time-dependent problems that require an initial condition 
for a specific solution. The (functional) density matrix is given by a superposition 
of wave functionals 



p(01,0 2 )=£pA(0l)*;(02) , (11) 
n 

where {^ n } is a complete set of wave functionals, and p n is the probability (J2 n Pn = 
1) that the system is in the state \l/ n . In general the ^/ n 's and the p n 's are time- 
dependent. 

In equilibrium the dynamics is time-translation invariant and energy is conserved. 
Then the complete set of wave functionals {^ n } can be chosen to be the set of the 
(time-dependent) energy eigenstates, while the p n 's are time-independent and are 
given by the canonical Boltzmann distribution, p n = e^'/fEn 6- ^")) where E n 
is the energy eigenvalue of the state \l/ n . The time evolution of the corresponding 
density matrix is trivial: it remains constant in time because the p n J s and the \l/ n \l/* 's 
are constant (N.B.: the time dependence of the \l/ n 's is just a phase). 

For non-equilibrium physics the time evolution of the density matrix is instead 
non-trivial. In fact, it might not be possible to choose the ^/ n 's as energy eigenstates, 
and the p n 's need not be Boltzmann factors and can change in time. Under the as- 
sumption that the time dependence of the \l/ n 's be determined by a time-dependent 
Schrodinger equation, the density matrix p satisfies the following differential equa- 
tion 

In order for Eq. ( |T2| ) to describe a well defined initial value problem for the time 
evolution of p, it is necessary to give the form of H and a model for dp n /dt. 

2 Related studies have been reported in Refs. jll], [l2|]; however, the formalism adopted there 
was somewhat different. In particular, the Liouville-vonNeumann equation (which in this section 
will arise only after a suitable approximation) was taken as a starting point. Moreover, as a 
result of the assumptions made about the geometry of the problem, previous DCC analyses using 
non-equilibrium quantum field theory are more relevant to the case of heavy-ion collisions. 



One important simplification is allowed when the evolution of interest is entropy- 
conserving; in fact, in this case one can work with time- independent p n 's, which 
indeed can be shown ||21|| to correspond to entropy-conserving time evolution. In 
these hypotheses, Eq. (|12|) takes the form of the quantum Liouville-vonNeumann 
equation with time-dependent Hamiltonian: 

While the investigation of the Liouville-von Neumann equation is still extremely 



difficult, it is certainly remarkably simpler than the original Eq. ([121). [We are 
not aware of any fruitful investigations of Eq. (|I2]).] It is therefore quite crucial to 
establish whether the assumption of entropy-conserving evolution (time-independent 
Pn's) is justified in the study of DCCs. In the baked Alaska scenario for the formation 
of DCC in hadron-hadron collisions that has been considered here, it appears that 
the assumption of entropy-conserving evolution could be quite accurate; in fact, 
the interior of the fireball is a rather "peaceful" place. Instead, in the case of 
heavy-ion collisions the relevant time scales are such that the formation of DCCs 
might be affected by non-entropy-conserving stages of the evolution of the quark- 
gluon plasma (i.e. one might be led to a definition of the system such that the 
evolution is entropy non-conserving). Even when it is safe to assume that the time 
evolution is entropy-conserving, and therefore the Liouville-von Neumann equation 
can be meaningfully taken as a starting point, for non-trivial Hamiltonians, such 
as the ones relevant for DCC physics, progress requires further approximations. 



Interestingly, one can obtain approximate solutions using the observation [21 



that the Liouville-von Neumann equation can be derived by varying an action-like 
quantity within a variational principle introduced by Balian and Veneroni [24}|. An 
approximate application of the variational principle, with a restricted variational 
ansatz, leads to approximate equations for the density matrix^. 



4 Outlook 

The recent increased interest in DCC physics has led to substantial progress, 
but much work is still needed. The challenge is quite strong both for theorists, 
who must find their way through the complicated models involved and provide 
experimentalists with a reasonable (quantitative!) picture of what to look for, and 
for experimentalists, who must disentangle the (possibly faint) DCC signal from 
substantial backgrounds. 

We are confident that our program of numerical and analytical investigation of 
the linear cx-model will soon lead to simulations realistic enough to be used in the 

3 Alternatively, as recently observed in Ref. (25) , an ansatz can be made directly at the level of the 
Liouville-von Neumann equation, without advocating the Balian- Veneroni variational riformulation 
of the problem. 



interpretation of hadron-hadron experimental results relevant for the DCC 
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